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OVERDENSITY OF /'-DROPOUT GALAXIES IN THE SUBARU DEEP FIELD: 
A CANDIDATE PROTOCLUSTER AT Z ~ 6 1 - 2 
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ABSTRACT 

We investigate the sky distribution of z ~ 6 Lyman break galaxies selected as i'-dropouts having 
i' — z' > 1.45 down to z' < 26.5 in the Subaru Deep Field (SDF). We discover 37 i'-dropouts clustered 
in a projected comoving 21.6 x 21.6 Mpc 2 region at z = 6, showing a local density excess. Carrying out 
follow-up spectroscopy we identify four of them as Lya emitters at z = 5.92, 6.01, 6.03 and 6.03 (spread 
over a distance of 46.6 Mpc). The number density of the cluster itself in SDF is ~ 2.2 x 10 -7 Mpc -3 , 
smaller than those of protoclusters (i.e., forming galaxy clusters) at z ~ 2-5.7. Also, the structure shows 
~ 4-21 times larger galaxy number density than those of z ~ 6 galaxies in a general field. It has a mass 
of M ~ 1.5J_ ; g x 10 15 M Q , comparable to those of z ~ 0-5 protoclusters. Since the contamination of our 
sample by interlopers is estimated to be quite low, 5.9%, most of the other unconfirmed i'-dropouts in 
the overdense region can be also z ~ 6 galaxies. Hence, it could be a candidate forming cluster at z ~ 6, 
representing a progenitor of galaxy clusters seen in the recent-day Universe. 

Subject headings: cosmology: observations — early universe — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

Observations of high-z star-forming galaxies such 
as Lya emitters (LAEs) and Lyman break galax- 
ies (LBGs) have been used to search for large- 
scale structures and primeval forming galaxy clusters — 
protoclusters — at high redshifts. Starting from the 
discovery of protoclusters mad e of these galax i es at 



2-3 dMalkan et all [19961 iPascarelle et all IJ 996: 



Giavalisco et all 11994. iLe Fevre et all 119961 : ISteidel et all 
1998), subsequent searches at higher redshifts have 



achieved remarkabl e success finding larg e -scale structures 
up to z ~ 4-5.7 (lOverzier et al.l 120061 : IVenemans et al 



2007t IShimasaku et al.l 120031: [Kashikawa et alj 1200 



% 



could be a protocluster at z ~ 6, and our follow-up spec- 
troscopy of some of the i'-dropouts forming this overdense 
structure. In §2, we describe the data used for our photom- 
etry as well as our z'-dropout sample selection. The sky 
distribution of z'-dropouts in SDF is analyzed statistically 
in §3. §4 presents the result of the follow-up spectroscopy. 
Finally, discussion and conclusion are given in §5 and 6. 
Throughout we adopt an (fl rn ,fl\,h) — (0.3,0.7,0.7) cos- 
mology and AB magnitude with 2" aperture unless other- 
wise specified. 

2. z'-DROPOUT SAMPLE 



Mal hotra et al.ll2005t lOuchi et al.l2005f) . Several attempts 



have been made to search for even higher- z galaxy clus- 
tering targeting fields around z ~ 6 q uasars, but only two 
possible such system has been f ound (S tiavelli et alj|2 005: 
IWillott et aLll2005t IZheng et alj|2006t lAiiki et al.ll2006f ). 

To study galaxy clustering at even larger lookback time 
(i.e., z > 6), we have selected z ~ 6 galaxy candidates by 
the z'-band dropout tec hnique, in the Subaru Deep Field 
dKashikawa et al.ll2004L SDF; 13 h 24 m 38. s 9, +27°29'25."9; 
an area of ~ 876 arcmin 2 ), taking the full advantage of 
its extremely wide and deep optical images. The present 
paper reports our finding of z'-dropout overdensity that 



2.1. Sample Selection 

We used official catalogs^ of the objects detected in 
SDF, for which each of the deep BVRi'z' images was 
taken by one pointing of the Subaru Prime Focus Camera 
(jMivazaki et alj|2002l Suprime-Cam) on the 8.2m Subaru 
telescope. All the images were convolved to have com- 
mon seeing size of 0."98. Source detection and photome- 
try in each wave band were performed b y using SExtrac- 
tor version 2.1.6 (|Bertin fc Arnoutall996f ) and the limiting 
magnitudes in our 2" aperture at 3ct are (B, V, R, i' , z') = 
(28.45,27.74,27.80,27.43,26.62). The 2" aperture and 
total magnitudes were measured with MAG_APER and 
MAG_AUTO parameters, respectively. From the ob- 

1 Based on data collected at Subaru Telescope, which is operated by National Astronomical Observatory of Japan. 

2 Some of the data presented herein were obtained at the W.M. Keck Observatory, which is operated as a scientific partnership among the 
California Institute of Technology, the University of California and the National Aeronautics and Space Administration. The Observatory was 
made possible by the generous financial support of the W.M. Keck Foundation. 

3 Cosmic Radiation Laboratory, RIKEN, 2-1 Hirosawa, Wako-shi, Saitama 351-0198, Japan 
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5 The Graduate University for Advanced Studies, 2-21-1 Osawa, Mitaka, Tokyo, 181-8588, Japan 
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jects detected down to z 1 < 26.5, we selected 158 i'- 
dropouts by using the criteria: B > 28.45, V > 27.74, 
i' - z' > 1.45. Null detections (< 3ct) in B and V 
were required since the flux of a LBG blueward of Lya 
should be absorbed by the intergalactic medium. Our 
i' — z' > 1.45 criterion picks up galaxies having Lyman 
breaks at z ~ 5-7.1, b ased on calculations using popu- 
lation synthesis models (jBruzual k. Chariot! 2003; M adaul 
Il995f h This model calculation also shows that there could 
be weak remainig fluxes not completely absorbed by the 
intergalactic medium in the i?-band wavelength of z ~ 6 
galaxies' spectra. Hence, we did not require a null detec- 
tion in R. All the i' magnitudes fainter than 28.62 (ler 
limits) were replaced by this value in application of the 
selection criteria. 

2.2. Contamination from Interlopers 

The problem in selecting z ~ 6 galaxies without near- 
infrared images is contamination. Three different types 
of interlopers could satisfy i' — z' > 1.45; dwarf star s 
with spectral type later than M7 (|Hawlev et al.l l2002f) . 
z ~ 1-2 old elliptical galaxies called Extremely Red 
Objects (EROs) with 4000A Balmer break in z'-band 
and z ~ 6 quasars with Lya absorption. We evalu- 
ated fractions of th ese interlopers i n our sample using 



the sa me method as Ota et all (|2005|) : using a star count 
model (|Nakaiima et al.l I2000T) . old elliptical ERO sam- 
ples actually observed in S ubaru/ XMM-Newton Deep Sur- 
vey (iMiyazaki et all 120031) and Hubble Ultra Deep Field 
(UP F) ( | Yan et aLll2004l), and a quasar luminosity func- 
tion (|Yan fc W indhorst 2004). The calculated number 
counts of each interloper as well as z ~ 6 LBG are shown 
in Figure Q] Most (~ 91%) of the i'-dropouts in SDF are 
fainter than z' = 25. The possible fractions of dwarfs, 
EROs and quasars in z 1 > 25 are only 5.9%, 0% and 
0.12%, respectiv ely. If we consider cosm ic variance of ERO 
numbe r density (ISomerville et al.1l2004[ ) and Poisson error 
for (jGehrelsl 119861) . the upper limit of ERO fraction 
could be < 16% in SDF survey area if th e se ER Os satisfy 
i' - z' > 1.45. However. [Malhotra et all (|2005l ) also con- 
firmed that no ERO shows i' — z' > 1.3 color (See their 
§3 and Fig.l) based on their spectroscopy of i-dropouts 
in UDF. Hence, we ignore the ERO contamination in the 
subsequent analyses and discussions. The low contamina- 
tion rate estimated here is consistent with the fact that 
our follow-up spectroscopy did not reveal any candidates 
to be interlopers (See @. 

3. SKY DISTRIBUTION 

3.1. Discovery of A Local Surface Overdensity 

The sky distribution of 158 i'-dropouts in SDF is shown 
in Figure [2j A very interesting region in the south shows 
a significant overdensity of these i'-dropouts which may 
be a protocluster at z ~ 6. To verify the significance 
of this compared to other regions quantitatively, we mea- 
sured the degree of density excess, a^2o, the overdensity of 
i'-dropouts in a circle of 20 Mpc comovin g radius at z = 6 
defined as of^ = [((£ - E) 2 } - E]/E 2 (| Peebles I [l980h . 
Here, E and E are the number of objects in the circle and 
the mean of E values measured in circles at 2000 randomly 
chosen positions in the z' image of SDF, respectively. We 
obtained a^20 = 0.63 with 3a significance for the region of 



our interest. This value is much higher than those of any 
other places in SDF, and comparable to that of z ~ 5.7 
protoclusters found bv lOuchi et al.l (|2005f ). This area also 
shows the highest density contrast: cts = (E — 'S)/'S = 3.6 
with 5cr excess measured in a comoving 8 Mpc circle. This 
value is comparable to those for a z ~ 4.8 structure found 
by Shimasaku et al. (2003, 2004) and z ~ 5.7 protoclusters 
bv lOuchi et alj (|2005f ). We found 56 and 22 i'-dropouts in 
20 Mpc and 8 Mpc circles, respectively. 

The overdense region contains two extremely bright 
iz' — 21.88 and 22.16) i'-dropouts that might be z ~ 6 
quasars or dwarf stars as seen in Figure [5] They could be 
part of a system consisting of LAEs physically associated 
with one or two quasars at z ~ 6, although these bright 
i'-dropouts are more plausibly dwarf stars. 

In the SDF, z ~ 6 galaxy candidat es had been also 
selecte d by two different techniques. Shimasaku et al.l 
(2005J) made zb and zr filters splitting the z' waveband 
at 9500A and selected out z = 5.9 ± 0.3 LBG candidates 
in the i' — zr vs. zb — zn colo r plane down to zr = 25.4. 
Meanwhile. fShiova et all J2005ft detected z = 5.8-6.5 LBG 
candidates (NB816-dropouts) using the i' — z' vs. NB816 
—z' color diagram down to z' = 26.07, where NB816 is a 
narrowband filter centered at 8150A with FWHM of 120A. 
Our overdensity region includes three such z^-detected ob- 
jects out of 14 i'-dropouts down to z' = 26.22 (total mag) 
assuming it is compara ble to zr = 25 A [z! ~ zr + 0.82 on 
average; see Table 2 in IShimasaku et all (|2005h ) and two 
NB816-dropouts out of 14 i'-dropouts down to z' — 26.07. 
They are regarded as more plausible galaxy candidates at 
z ~ 6. Note that since their candidate selections adopt 
null detection in i?-band (R < 2-3cr), we also limited our 
sample to R < 2-3ct in the comparison. 

4. FOLLOW-UP SPECTROSCOPY 

To obtain information about the physical clustering and 
to check the validity of our sample selection, we obtained 
spectra of 20 i'-dropouts in the overd ense region with 
Faint Object Camera and Spectrograph (jKashikawa et all 
12001 FOCAS ] on the Subaru T elescope and the DEIMOS 
spectrograph (jFaber et al.ll2003l ) on the Keck II Telescope. 
Spectrosc opy of these ^-dropouts is summarized in Table 
H and in iKashikawa et alj (J2006D . live et all (|2006j) and 
Otaet alJ|2008l) . 



We identified four i'-dropouts as LAEs at similar red- 
shifts z = 5.92, 6.01, 6.03 and 6.03. The emission lines 
found in their spectra are asymmetric with their blue 
side absorbed and have weighted skewness of Sw = 3.43- 
17.9A, comparable to the values of z — 5.7 and 6.6 LAEs, 
S^ 5 - 7 = 3 33-20.66A and ffT 6 - 6 = 3.12 -16.34A (See Ta- 
ble [21 and IShimasaku et al.1 (|2006f ) and IKashikawa et al.l 
(2006)). Since each of the four spectra shows only one 
emission line, the line is not H/3, [OIII]4959A, [OIII]5007A, 
Ha, [SII]6717A or [SII]673lA. Also, the emission lines we 
found are not the z — 1.6 [Oil] doublet with 7A separation 
(3726A and 3727A at rest frame) because we actually re- 
solved even closer (5-6A) skylines in the same wavelength 
range. In addition, the objects are not z ~ 6 active galac- 
tic nuclei, since we did not find Ny 1240A lines at 8680A. 
Also, FWHMs of the 4 confirmed emssion lines are 367- 
464 km s _1 and much smaller than the avera ge FWHM 
of z ~ 6 quasars' Lya lines, ~ 6,000 km s _1 (JFan et all 
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12004ft . Hence, we consider all the emission lines we de- 
tected to be Lyo of LAEs. The photometric and spectro- 
scopic properties of four confirmed LAEs as well as their 
spectra, are presented in Table [2] and Figure [3] Their sky 
distribution is shown in Figure O 

Spectra of 15 out of 16 unconfirmed i'-dropouts were ob- 
tained during the 2005 run whose main targets were other 
brighter objects, and thus integration times were rather 
short (only 1.25-3.75 hours) for our faint (z' = 25.5- 
26.5) i'-dropouts. Hence, if the Lya lines are absent or 
overlapped with OH-airglow lines, identification is much 
harder. Thus, we cannot completely deny the possibility 
of these unconfirmed i'-dropouts being galaxies at z ~ 6. 
None of the 20 spectroscopic targets turned out to be low- 
z emission-line galaxies. Hence, the possible contamina- 
tion of low-z galaxies in our photometric i'-dropout sample 
seems very rare. 

5. DISCUSSION 

Analyses of the overdense region implies that it might 
be a z ~ 6 protocluster. In the following discussion, we 
estimate the basic properties of the overdense region and 
compare them with those of z < 6 protoclusters and field 
galaxies to see if it can be a z ~ 6 protocluster that is a 
progenitor of galaxy clusters seen in the present-day uni- 
verse. 

5.1. Number Density of the Cluster Itself 

Since we have probed the comoving depth of 887 Mpc 
(corresponding to z ~ 5-7.1 of z'-dropout selection func- 
tion) and area of SDF (876 arcmin 2 ), our survey volume 
is approximately 4.5 x 10 6 Mpc 3 . Since we found only 
one overdensity in this volume, its number density is es- 
timated to be 2.2 x 10~ 7 Mpc -3 . This is smaller than 
those of protoclusters at lower-z; e.g., 2.2 x 10~ 6 Mpc~ 3 at 
z = 5.7 dOuchi et al.ll2005ft . 7.1 x 10~ 6 Mpc" 3 at z = 4.86 
(Shimas akuet alJl2003h ,~~ 6.0 x 10~ 6 Mpc" 3 at z = 2-5 
( Venemans et al.M2007f) . ~ 3.0 x 10~ 6 Mpc" 3 at z ~ 3 
( Stci del et al.l Il998t converted to our cosmology) . The 



smaller number density might be because forming clusters 
are not virialized yet and difficult to identify at higher 
redshifts as implied by hierarchical clustering scenario. 

5.2. Galaxy Number Density 

If we define the border of the overdense region as the 
place where the density contrast drops to the half of unity 
5s ~ 0.5, its size is ~ 9' x 9' (21.6 x 21.6 comoving 
Mpc 2 at z = 6), approximately within ARA = 13-22' 
and ADEC = l'-10' in Figure |2J There are a total of 37 
i'-dropouts in it. Four of them are identified as LAEs, of 
which one has z' < 25. Also, two bright objects might be 
dwarf stars or quasars. Since the contamination rate of 
dwarfs in z' > 25 is expected to be 5.9% (See j g~IJ) , 32 out 
of the remaining 34 i'-dropouts in z' > 25 can be z ~ 6 
galaxies. Hence, we assume that the overdensity consists 
of 33 galaxies at z ~ 6. If we assume the line-of-sight 
scale of the overdensity corresponds to Z = 5.92-6.03 of 
confirmed LAEs (46.6 comoving Mpc), its volume would 
be 2.2 x 10 4 Mpc 3 . This gives the galaxy number density 
of "cluster ~ l-S^Lo x 10~ 3 Mpc -3 . The uncertainty is 
large since it includes Poisson errors as well as large cosmic 



variance due to the small volu me. Cosm i c vari ance is esti- 
mated in the same way as in lOta et all (|2008ft . For com- 
parison, we calculated the number density of z ~ 6 galaxies 
in the general field and obtained nfidd ~ 4.0^4 x 10 -4 
Mpc -3 (same uncertainty as above) by inte grating the 
best- fit Schechter luminosity fun ction (LF) (ISchechter I 
Il976ft of i-dropouts derived by Bou wens et al.l (J2006) down 
to Mi5oo,AB = —20.2 (corresponding to our detection 
limit z' = 26.5). Thus, the overdensity region has about 
four times (riciuster/ftfield ~ 3.8^ 2 'e) higher surface density 
than galaxies in a field. Note that this factor and n c i us ter 
are lower limits because the incompleteness regarding i'- 
dropout selection is not taken into account. Hence, those 
values could be larger if the incompleteness were corrected. 

5.3. Star Formation Rate 

We estimate the average UV star formation rate (SFR) 
of the overdense region, excluding two bright i'-dropouts 
that might be dwarfs or quasars and assuming that all the 
rest of 33 z'-dropouts (corrected for contamination) are 
galaxies at Z = 6. Their rest frame UV continuum lumi- 
nosities, L„, are calculated from z'-band total magnitudes 
and converted to S FR at 1500A, using the formalism of 
iMadau et al.1 (| 1998ft : L v = 8.0 x 10 27 SFR erg s" 1 Hz" 1 . 
This yields the average SFR of 11.5-M0yr — x . Meanwhile, 
the average SFR of 116 z'-dropouts with the same magni- 
tude cuttoff of z' > 24.5 (total mag) outside the overdense 
region in SDF is 12.0M Q yr _1 . This implies that star for- 
mation activity in the the overdense region does not differ 
much from that in a general field. 

5.4. Mass Estimate 

Now that the galaxy overdensity 5 ga \ = n c i us tcr/^ncid — 
1 ~ 2.8^26 i s known, we can estimate the mass of the 
overdensity region with M = pV(l + 8 m ) (jSteidel et al.l 
ll998HVenemans et al .112005ft . Here, p = n m (3H /8irG) 2 = 
4.1 x 1O 1O M0 Mpc -3 is the mean matter density of the 
Universe with gravitational constant G. V = 2.2 x 10 4 
Mpc 3 is the volume of the candidate protocluster. 8 m is 
the mass overdensity related to the bias parameter b with 

1 + bS m = C(l + <5 gal ). C = [1 + / - /(l + SU 3 )}, where 
/ ~ 1 in our cosmology, corrects the volume for redshift 
space di stortion due to peculiar motions of prot ocluster 
galaxies (jSteidel et al.lll998l : IVenemans et al.ll2005ft. Using 
<5 ga i ~ 2. 8^2:6 and & = 3.4 ± 1.8 suggested by lOuchi et al.l 
(|2005ft . we estimate that S m - 0.62tgg and thus M ~ 
l.SJIo I x 1O 15 M0. Our estimated mass is comparable to 
the mass of the z = 4.1 protoclust er, (1-2) x 10 15 M Q , 
estimated by IVenemans et al.l (J2002ft , that of the z ~ 3 
LBG c oncentration, (0.98-1.4) x 1O 15 M bv lSteidel et~aT1 
(11998ft. those of z = 2-5 protoclusters, (2-9) x 10 14 M Q , 
by IVenemans et all (|2007l ), and that of a z = 4.8 6 pro- 
tocluster, 3 x 10 14 M Q , bv IShimasaku et al.l (|2003ft . Ac- 
cording to the large sample statistics (152-1 70 clusters) of 
nearby jz < 0.15) galaxy clusters studied bv lGirardi et al.l 
( 1998a. b), our overdensity mass is also comparable to their 
masses. Moreover, the estimated mass and number den- 
sity of the overdensity measured in M5.1l are consistent with 
mass functions derived from large samples of ~ 72-300 
z ~ 0.1-0.2 nearby galaxy clusters l|Girardi et al.l [l998b: 
Bah call et al. 2003; Rines et al.l l2007f ). Hence, like other 
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lower- z protoclusters, it might be an ancestor of a cluster 
of galaxies seen in the present-day Universe. 

6. CONCLUSION 

The number density, number density of member galax- 
ies and mass of the i'-dropout overdensity in SDF imply 
that there exist forming clusters — protoclusters — at a red- 
shift as high as z ~ 6. However, the result presented here 
is based on photometry and a few spectroscopically iden- 
tified LAEs at z = 5.92-6.03. Hence, the overdensity is 
still a "candidate" protocluster. Further follow-up deep 
spectroscopy of the unconfirmed i'-dropouts will reveal if 



it is a real z 
properties. 



6 protocluster and constrain its physical 
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Table 1 
Status of the follow-up spectroscopy of i'-dropouts in the overdense region 
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grating 


R 


wavelength 


filter 


slit 


standard stars e 






(HST) 


(") 


(seconds) 




(mm - 1 ) 




(nm) 




(") 




2 a 


1 


24-27 Apr 2004 c 


0.4-0.8 


16200 


F 


300 


1300 


540-1000 


058 


0.6 


Hz 44, Feige 34 


8° 





14-15 May, 1 Jun 2005 d 


0.5-1.0 


13500 


F 


175 


1600 


830-1000 


SDSS z' 


o.s 


Hz 44, Fiege 34, Feige 110 


1 


1 


14-15 May, 1 Jun 2005 d 


0.5-1.0 


9000 


F 


175 


1600 


830-1000 


SDSS z' 


o.s 


Hz 44, Fiege 34, Feige 110 


5 





14-15 May 2005 d 


0.5-1.0 


4500 


F 


175 


1600 


830-1000 


SDSS z' 


o.s 


Fiege 34, Feige 110 


2 


2 b 


23-24 Apr 2004 c 


0.6-1.0 


9000 


D 


830 


3600 


500-1000 


GG495 


1.0 


Feige 110, BD+28°4211 



NOTE: Columns in the table show the number of i'-dropouts spectroscopically observed (Ob), the number of i'-dropouts identified (ID), 
the dates of observation, typical seeing size during the observation, exposure time, the instrument (I) used (F for Subaru/FOCAS and D for 
Keck/DEIMOS, respectively. 

a One of them was observed in both dates. 

b One of them was identified bv lNagao et al.l (|2005l . I2007t l. 

c Spectra were obtained in the observations by Kashikawa ct al. (2006) 

d Spectra were obtained in the observations by live et alj l|2006f) and lOta et al.l (2008). 

e See lOkel fi990T) and lHamuv et al.l Jl99^) for details. 
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Table 2 
Properties of four confirmed z ~ 6 Lyq emitters in the overdense region 



Object and Coordinate 


z 


i> 


z> 


S m 


F(Lya) 


L(Lya) 


SFR(Lya) 


FWHM 


EW (Lya) 






(mag) 


(mag) 


(A) ( 


10~ 17 crg s" 1 cm -2 ) 


(lO^ergs- 1 ) 


(M s yr" 1 ) 


(A) 


(A) 


(1) SDF J132418.4+271633" 


5.92 


27.30 


25.74 


3.43 ± 1.07 


1.3 + 0.2 


0.49 + 0.07 


4.43 + 0.6 


13.5 


36 


(2) SDF J132429. 0+271918° 


6.01 


>28.62 


26.29 


8.52 ± 1.39 


0.84 + 0.5 


0.34 + 0.19 


3.06+ 1.7 


11.2 


30 


(3) SDF J132426.5+271600° ,c 


6.03 


27.43 


25.36 


5.12 ± 1.47 


3.6 + 0.3 


1.5 + 0.1 


13.7 + 0.9 


13.2 


84 


(4) SDF ,1132431.5+271509° 


6.03 


27.83 


25.88 


17.9 + 0.32 


4.3 + 0.3 


1.7 + 0.1 


15.6 + 0.9 


14.2 


81 



NOTE: Units of coordinate are hours: minutes: seconds (right ascension) and degrees: arcminutes: arcseconds (declination) using J2000.0 
equinox, i' and z' are 2" aperture magnitudes. Magnitude is replaced by its l a limit if it is fa i nter t han the limit. The indicato r of a 
line profile asymmetry, weighted skewness S w , is calculated using the definition in Shimasaku ct al. (2006) and Kashikawa ct al. (2006). The 
Lyo line fluxes, i^Lya) measured in the spectra are used to calculate the Lyct line luminosities, L(Lya). Then, L(Lya) is converted into 
the corre sponding star for mation rates, SFR(Lya), using the relation, S'.F^Lya:) = 9.1 x 10 _43 L(Lya)AfQyr — 1 , derived from Kennicutt's 
equation (Kennicutt 1998) with the case B recombination theory. The FWHM of each Lyo li ne is measured by fi tting a Gaussian profile. The 
rest frame equivalent widths of Lyo lines EWo (Lyo) are calculated in the same manner as in INagao et al.l l)2007l) . 

"Confirmed by Subaru/FOCAS. 

"Confirmed by Keck/DEIMOS. 

identified bv lNaeao et~aT1 [120051 , 120071) . 
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Fig. 1. — Differential number counts of all the possible i'-dropout objects with z' > 23 (total magnitude) in SDF. They are shown in 
several different symbols and lines as labeled in the figure. The vertical dotted line is the 3<r limiting magnitude of z' = 26.62. All the vertical 
errors are Poissonian errors. Poissonian errors for dwarfs and total i'-dropout counts are omitted for clarity. The count of z ~ 6 LBGs is 
the difference between i'-dropouts total and the sum of M/L/T dwarf star and z ~ 6 quasar. ~ 91% of the i'-dropouts have z' > 25 where 
the contamination rates of dwarfs, EROs and quasars are only 5.9%, 0% and 0.12%, respectively. The surface density of i'-dropouts in the 
9' X 9' overdense region (See N5.1l for the size definition) is also shown by asterisks. It is a factor of 2—3 times larger than the surface density 
of i'-dropouts in the entire SDF field. 
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Fig. 2. — Sky distribution of i'-dropouts in the Subaru Deep Field. They are filled circles with sizes proportional to their z'-band total 
magnitudes (huge red: z' < 24.5, very large orange: 24.5 < z' < 25, large magenta: 25 < z' < 25.5, medium green: 25.5 < z' < 26, small 
blue: 26 < z < 26.5). North is up, and east is to the left. The comoving scale projected to z = 6 is also shown along the axes. Solid lines 
are density contours of &■£ = 0, 0.2, 0.4, 0.6, 0.8 and 1.0 (See §3.1 for the definition of 5^). Four spectroscopically confirmed Lyct emitters 
at z = 5.92, 6.01, 6.03 and 6.03 are encircled and labeled with numbers 1—4 as in Table [2] and Figure [3] They are located in the overdense 
region to the south. Two i'-dropouts enclosed by squares are extremely bright (z' = 21.88 and 22.16) and might be z ~ 6 quasars physically 
associated with z ~ 6 galaxies or possibly dwarf stars. Three and two i'-dro pouts enclosed by diamo nds and triangles in the overdense region 
are more plausible c andidates fo r z ~ 6 LBGs selected by using zr filter by Shimasaku et al. (2005J) and by the narrowband NB816-dropout 
technique by Shi ova et al.l 1120051) . 
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Fig. 3. — Spectra of four confirmed z ~ 6 Lya emitters in the overdense region as well as OH skylines taken with Subaru/FOCAS in the 
same wavelength region. Spectra obtained using Subaru/FOCAS and Keck/DEIMOS are labeled (F) and (D), respectively. The labels of 
(l)-(4) correspond to those in Figure[2]and Table|2] 



